function RX = Channel(fd,f0,TX,SampleOffset,t,snr,fc)
    Doppler_phase=exp(1j*2*pi*fd*t);
    Fixed_phase=exp(1j*2*pi*f0*t);
    Total_phase=exp(1j*2*pi*(f0+fd)*t);
    %% 双径
    TX1=TX.*Total_phase;
    TX2=TX.*Fixed_phase;
    RX1=[zeros(1,SampleOffset),TX1];
    RX2=[zeros(1,floor(SampleOffset/2)),TX2];
    Dist_Real=100;
    Dist_BS=50;
    PathLoss1=32.44+20*log10(Dist_Real/1000)+20*log10(fc/1e6);
    PathLoss2=32.44+20*log10(Dist_BS/1000)+20*log10(fc/1e6);
    RX1_noise=awgn(10^(-PathLoss1/10)*RX1,snr,'measured');
    RX2_noise=awgn(10^(-PathLoss2/10)*RX2,snr,'measured');
    RX=awgn(RX1(1:length(RX2))+RX2,snr,'measured');
end